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ABSTRACT 



The G9.5 giant e Oph shows evidence of radial p-mode pulsations in both radial velocity and luminosity. We re-examine the observed 
frequencies in the photometry and radial velocities and find a best model fit to 18 of the 21 most significant photometric frequencies. 
The observed frequencies are matched to both radial and nonradial modes in the best model fit. The small scatter of the frequencies 
about the model predicted frequencies indicate that the average lifetimes of the modes could be as long as 10-20 d. The best fit model 
itself, constrained only by the observed frequencies, lies within ±lcr of e Oph's position in the HR-diagram and the interferometrically 
determined radius. 
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1. Introduction 

In radial velocity measurements from ground, De Ridder et al. 
(2006) found power excess in the frequency spectrum of e Oph 
with a maximum amplitude of about 3.5 ms"' at a frequency of 
approximately 60/iHz (~5.2 d '). However, the modest duty cy- 
cle of the groundbased observations and the resulting d"' aliases 
made it difficult for them to choose between possible p-mode 
spacings of ~5 and ~1 //Hz. Barban et al. (2007) analyzed re- 
cent MOST (Microvariability & Oscillations of STars) satellite 
observations and obtained frequencies, amplitudes, and mode 
lifetimes for a sequence of radial modes. Their mode lifetimes 
(~2.7 d) compared to another giant ^ Hya are consistent with the 
lifetimes determined from observations (Stello et al. 2006) but 
are shorter than the lifetimes inferred from theoretical consider- 
ations (e.g., Houdek and Gough 2002). 

Barban et al. (2007) only considered radial modes in their 
analysis. Indeed, it has been suggested that only radial /5-modes 
should be observable in giants owing in part to the complex 
frequency spacings of nonradial mixed (or bumped) modes. 
The theoretical position is unclear. The nonadiabatic models 
of Dziembowski et al. (2001) show the existence of strongly 
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trapped unstable (STU) nonradial modes in the central regions of 
giants in between very closely spaced modes of the same degree 
/. The STU modes follow the same pattern of frequency spac- 
ing as radial acoustic modes (i.e., p-xnodes). They correspond to 
mixed ^-modes//?-modes. Furthermore, the authors find more / 
- 2 STU modes to be unstable than I = 1 STU modes for stars in 
the luminosity range of e Oph (see Fig. 4 in Dziembowski et al. 
2001). Only in their linear stability analysis, which includes the 
effects of turbulent pressure based on the mixing-length formal- 
ism, do they find predicted amplitudes of the radial modes of the 
sub-giant star a UMa to be much smaller than observed, with the 
nonradial mode amplitudes predicted to be even smaller. On the 
observational side, Hekker et al. (2006) note that the line pro- 
file variations of several pulsating red giants (including e Oph) 
suggest the existence of nonradial modes. 



If some or all pulsating red giants do have radial and observ- 
able nonradial modes with relatively long lifetimes, then these 
modes will enable asteroseismologists to constrain the deep in- 
teriors of red giant stars, as well as set limits on the excita- 
tion mechanisms themselves. We are therefore motivated to re- 
examine the MOST photometry of e Oph, combined with the ra- 
dial velocities obtained by De Ridder et al. (2006), to explore the 
possibility of longer-lived (longer than 3 days) nonradial pulsa- 
tions. 
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Fig. 1. MOST photometry of e Oph (V=3'"24). For better visibil- 
ity, circles are averages of 20 consecutive data points (covering 
about 10 min each) of the smoothed raw data (running average 
of 350 exposures); the un-averaged measurements of e Oph are 
shown as small grey dots. The inset is an enlargement of 4 day 
subset of the data. 



2. Observations 

2.1. Spectroscopy 

The radial velocity data of e Oph were provided to us by De 
Ridder et al. (2006) from their double-site campaign with the 
fibre-fed echelle spectrographs CORALIE (Swiss 1.2-m Euler 
telescope at La Silla, Chile) and ELODIE (French 1.93-m tele- 
scope at the Observatoire de Haute-Provence). The spectra were 
obtained during 53 nights (spanning 75 days) in summer 2003, 
with S/N of at least 100 in every spectrum at a wavelength of 
550 nm. Radial velocities were measured from the spectra us- 
ing the optimum-weight method relative to a nightly reference 
spectrum. This suppresses any signal (including signal intrinsic 
to the star) with periods longer than about 1 day. 

2.2. Photometry 

MOST houses a CCD photometer fed by a 15-cm Maksutov tele- 
scope through a custom broadband optical filter (350 - 750nm). 
The satellite's Sun-synchronous 820km polar orbit (period - 
10 1.4 min) enables uninterrupted observations of stars in its 
Continuous Viewing Zone (+36°>5>-18°) for up to 8 weeks. 
A pre-launch summary of the mission is given by Walker et 
al. (2003) and on-orbit science operations are described by 
Matthews et al. (2004). 

The MOST photometry of e Oph is described by Barban et 
al. (2007). The data cover 28.4 d during 16 May - 13 June 2005, 
and were obtained in Fabry Imaging mode, where the telescope 
pupil, illuminated by the star, is projected onto the CCD by a 
Fabry microlens as an annulus covering about 1300 pixels. To 
meet satellite downlink restrictions, the data are downloaded in 
two formats: resolved CCD subrasters of the stellar and back- 
ground Fabry images (Science Data Stream 2 or SDS2), and 
compressed data preprocessed on board the satellite (SDS 1). We 
use the resolved SDS2 data for our analysis here because it al- 
lows us to carry out background decorrelation of the Fabry im- 
ages (Reegen et al. 2006). Our subset of the MOST photometry 



contains 83,238 measurements with an exposure time of 7.5s, 
sampled every 30s. 

As in the Barban et al. (2007) analysis, for our re-analysis of 
e Oph we reject images obtained during MOST orbital phases 
of high scattered Earthshine, which for the e Oph run, occurred 
when the satellite was passing over the fully illuminated North 
Pole. We introduce gaps in the data of about 24 min duration 
every 101. 4 min (the MOST orbital period). We also reject all 
3cr outlier frames. In total 24% of the exposures were rejected 
yielding an overall duty cycle of 76%. The regular gaps per orbit 
introduce aliases in the spectral window of the data, but these 
alias peaks are spaced ~14.2d"' (~165 fiHz) apart. Fortunately, 
this alias spacing, unlike the Id"' ground based alias spacing, 
does not introduce any ambiguities in the frequency identifica- 
tions in the range where e Oph exhibits p-mode pulsations. 

MOST was not originally designed to obtain differential pho- 
tometry, but thanks to its stable thermal environment and com- 
pact design, the instrument is sensitive to stellar variability at the 
millimagnitude level (and below) over timescales of up to one 
week. Only for very long timescales, comparable to the duration 
of the observations, instrumental drifts became severe and super- 
pose possible intrinsic variations. Stars like the red giant e Oph 
have turbulent envelopes which probably introduce longer-term 
variability than the /7-modes. Regardless, all long term variations 
(i.e., longer than a week), whether or not they are intrinsic, are 
ignored. 

Fig. [T] presents the binned light curve of e Oph, as well as 
the original unbinned reduced photometry. The data show very 
complex variations on timescales of hours to days, also seen by 
Barban et al. (2007), in addition to some slower trends. The clean 
spectral window of the data, with alias peaks spaced by 14. 2d"', 
shows that the slower trends do not interfere with the /5-mode 
identifications. There exists some modulation of stray light at 
low harmonics of Id"' . This is caused by the fact that the MOST 
orbit carries the spacecraft over nearly the same location on the 
Earth (hence, a similar albedo feature) every day. 

The reduced MOST e Oph photometry used in this paper, as 
well as the raw data and the Barban et al. (2006) reduction, are 
available in the MOST Public Data Archive on the Science page 
at www.astro.ubc.ca/MOST. 

3. Frequency Analysis 

The Fourier amplitude spectra of the e Oph MOST photometry 
(already prewhitened for instrumental signal; see latter in this 
section) and the radial velocity (RV) data are plotted in Fig.|2] 
along with the spectral window functions of both data sets. In 
the photometric amplitude spectrum, there are no aliases in the 
frequency range plotted, since the sampling sidelobes fall well 
outside the p-mode frequency range. There are, however, peaks 
at 2 and 3d"' (~23.1 and ~34.7yuHz) due to modulation of stray 
light in the data, as mentioned above. Some of the power at fre- 
quencies below 10 fjHz may be due to intrinsic turbulence in 
the stellar envelope, and some may be due to slow instrumental 
drifts but there is no significant leakage of this power to higher 
frequencies in the p-mode range because of the clean spectral 
window. 

The RV amplitude spectrum shows the cycle/day (d"') alias 
patterns due to the daily gaps in the groundbased spectroscopy. 
The data were subjected to a high-pass filter by De Ridder et 
al. (2006) so amplitudes at frequencies below ~20//Hz are sup- 
pressed. The RV spectrum has better frequency resolution than 
the photometric spectrum because of the longer time coverage 
of the RV data. 
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Fig. 2. Top panel: Fourier amplitude spectra of e Oph, and that of an intrinsically constant MOST guide star, HD 146490 (inverted 
for easier comparison). Both are shown at the same scale. Middle and bottom panel: Residual amplitude spectra of e Oph MOST 
photometry after prewhitening instrumental signal and radial velocity data (De Ridder et al. 2006). The spectral window functions 
of both data sets are shown in the insets; note the truncated amplitude scale for MOST spectral window. 



We used the routine SigSpec (Reegen 2007) to identify sig- 
nificant frequencies in the Discrete Fourier Transform (DFT) 
spectra of both data sets. SigSpec employs an exact analytical 
solution for the probability that a peak of given amplitude could 
be generated by white noise. Its main advantage over commonly 
used signal-to-noise ratio estimates is that SigSpec appropriately 
incorporates the frequency and phase angle in Fourier space and 
the time domain sampling. On average, a signal-to-noise ratio of 
4, suggested as a reliable significance estimator by Breger et al. 
(1993), roughly corresponds to a SigSpec significance value of 
5.46. The SigSpec significance is equal to logic the probability 
that the peak is not due to white noise. In other words, a signal 
amplitude of 4 times the noise level would appear by chance at 
its given frequency in only one of 10^^*" cases, assuming white 
noise. 

SigSpec automatically identifies the peak with the largest sig- 
nificance in the frequency range searched, determines the am- 
plitude and phase associated with this frequency through least- 
squares sinusoidal fitting, and subtracts that signal from the time 
series. The residuals are then used in the next SigSpec signif- 
icance calculation, and the process is repeated until there are 
no frequencies remaining with SigSpec significance above the 
threshold value set for the search. 

In the MOST photometry we find 99 peaks with SigSpec sig- 
nificances (sig) greater than 4.8 (equivalent to S/N ~ 3.75) from 
- 350;uHz (~30d '). 17 of the peaks are coincident with the 
MOST satellite orbital frequency (~165/iHz or 14. 2d"') and its 



first overtone as well as with sidelobes produced by the 1 cy- 
cle/day and long term modulation of the stray light amplitude. 
We compare the remaining 82 frequencies (all below ~120/iHz 
or ~10d"') with frequencies in the background signal and with 
frequencies in the light curve of the simultaniously observed 
guide star HD 146490 (V = 7^2, sp. A2). (MOST guide star 
photometry is collected from defocused star images in subrasters 
of 40x40 pixels. Although, not suitable for diflFerential compar- 
isons with the Fabry Imaging photometry, the guide star photom- 
etry can be used to help identify instrumental or environmental 
frequencies in the data.) The background time series was ex- 
tracted from the raw e Oph Fabry Image data using the average 
intensity of pixels well outside the annulus illuminated by stel- 
lar light that contains only signal from the sky background along 
with the bias and dark current levels of the CCD. SigSpec iden- 
tified 34 peaks, from to 120;uHz, with sig > 4.8 in the back- 
ground signal, and 14 peaks in the guide star, HD 146490, time 
series. 

We reject all frequencies that occur (within the uncertainties) 
in both the e Oph photometry and either the background or HD 
146490 time series, even if the amplitudes in the latter data sets 
are small. The frequency uncertainties are estimated according 
to Kallinger et al. (2007b) who define an upper limit for the fre- 
quency error of a least-squares sinusoidal fit as (T ■ V^) ' for 
a monoperiodic signal and (4T)"' if a second signal component 
is present within twice the classical frequency resolution defined 
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by the time series length T. We are left with 59 frequencies rang- 
ing from 3.1 lyuHz to 93.7yL(Hz (see Tab.[Tli. 

Figure|2] illustrates our problem to safely determine intrin- 
sic frequencies below approximately 20/iHz. The average am- 
phtudes increase by a factor of about 3 and it is not clear - 
nor can it be decided on the grounds of the present photom- 
etry - which contribution comes from possibly yet undetected 
instrumental properties, granulation signal (1/f characteristics), 
pulsation, artefacts enforced by short mode lifetimes of genuine 
signal, or a combination of all this. Furthermore, SigSpec - de- 
signed to handle white noise only - cannot be applied anymore 
in this low frequency region. An unsatisfying work-around to 
take the increasing noise level towards the low frequencies into 
account could be to increase the acceptance level in the low fre- 
quency region from sig = 4 to about 36, as sig scales with power. 
We divide Tab.[T]in a set of frequencies below 20//Hz, not used 
for modeling e Oph, and above 20//Hz. 

We are then left with 39 frequencies (20 to 93.7/vHz). 
Because it is extraordinarily difficult to fit a model to a set of 
frequencies when the list of frequencies contains even a small 
number of spurious frequencies, i.e., of non stellar origin, we re- 
strict our analysis to only those frequencies above 20juHz with 
sig > 10.0 (roughly corresponding to amplitudes > 60 ppm). We 
are convinced that these large amplitude frequencies are truly 
intrinsic and not contaminated by instrumental or other effects 
(different from pulsation). There are 21 frequencies that satisfy 
these conditions, more than enough to constrain the models. The 
low frequency set was not used for modelling, but is listed for 
completeness. However, a surprinsingly large fraction of fre- 
quencies (8 of 19, indicated by an asterisk in Tab.[T]l coincide 
with predictions of our best model (see Sec. 4.1). 

We searched the RV data for frequencies from - 120yuHz 
(~10d"') down to a sig limit of 5.2 (corresponding to S/N ~ 3.9). 
We identify 25 frequencies. These are more ambiguous identifi- 
cations because of the aliasing of the RV data; some of the fre- 
quencies may be ±1 d"' sidelobes of the real values. 11 of the 
25 frequencies identified in the RV data are also present, within 
the frequency uncertainties, in the MOST photometric data. All 
significant RV frequencies are listed in Tab.[T] 

4. Asteroseismic Analysis 

As a first step toward understanding the nature of the oscilla- 
tions in e Oph, we attempt to match the observed frequency 
spectrum with the eigenspectra determined from stellar mod- 
els. We calculated radial and nonradial pulsation frequencies 
of models distributed along evolutionary tracks on the giant 
branch of the HR Diagram. The models were constructed with 
the Yale Rotating Evolutionary Code (YREC, Guenther et al. 
1992), adopting the OPAL equation of state (Iglesias & Rogers 
1996), OPAL opacities in the interior (Rogers & Iglesias 1994) 
and the low-temperature opacities from Alexander & Ferguson 
(1994) in the envelope. Diffusion of helium and heavy elements 
was not included. The model eigenspectra were generated by 
the nonradial, nonadiabatic stellar pulsation program developed 
by Guenther (1994). The program solves the linearised nonadia- 
batic oscillation equations using the Henyey relaxation method. 
The inert helium core in giant stars supports closely spaced 
g-modes at frequencies that overlap the frequency range of the p- 
modes. As a consequence there is significant mode mixing (also 
known as mode bumping) where the nonradial /?-modes and g- 
modes interfere with each other (see Kallinger et al. 2007a and 
Guenther et al. 2000). Any nonradial /?-modes seen on the sur- 
face of the star have g-mode character in the interior. For gi- 
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Fig. 3. Echelle diagram of the 21 most significant observed os- 
cillation frequencies (above 20yuHz) used to constrain the fit and 
the best fitting model modes. The observed photometric modes 
are shown as circles with error bars indicating observational un- 
certainties. Filled circles show to the 18 observed frequencies 
matched by model modes and open circles to the remaining 3 
frequencies. The model modes are given by grey symbols con- 
nected by line segments. 



ants the radial order n of the g-mode part of the mode can be 
very high (~20 to 100). These are the same modes identified as 
strongly trapped unstable (STU) modes by Dziembowski et al. 
(2001). The densely packed fluctuations in the radial displace- 
ment eigenfunctions near the stellar core heavily strain the nu- 
merical calculation and can lead to numerical instabilities that 
prevent convergence of a solution. To avoid this complication, 
we force the code to ignore the inner 10% (in radius) of the stel- 
lar models. We have shown (Kallinger et al. 2007a) that the mode 
eigenfunctions in the outer 90% of the model are practically the 
same if the inner part of the model is used or not. The disad- 
vantage of this method is that one does not obtain the correct 
number of radial nodes for the ^-mode part of the eigenfunction. 
For more details, see Kallinger et al. (2007a). 

We computed radial and nonradial eigenfrequencies (with 
mode degree / up to 3) for more than 37,000 models of 
intermediate-mass stars (from 1.5 to 3.5Mq in steps of 0.02Mq) 
during their H-shell burning phase. The model sequences start 
at temperatures around 5500K (near spectral type Gl) and con- 
tinue up to the initiation of He-core burning at the end of the red 
giant phase. We set the initial hydrogen and metal mass fractions 
to (X, Z) = (0.71, 0.01), to match the observationally determined 
metal abundance from Cayrel de Strobel et al. (2001), who found 
[Fe/H] = -0.25. The mixing length parameter, which adjusts the 
temperature gradient in convective regions, was set to a = 1 .74 
to match the findings of Yi et al. (2003). For the model physics 
we are using, a calibrated solar model would have a - 1.63. The 
value of a we adopt yields a model at the solar age and luminos- 
ity that is hotter than the Sun by ~30K. 
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The best fits between the observed frequency spectrum and 
our grid of model frequencies were found by searching for min- 
ima in;^'^, defined as: 



X' 



,2 _ J_ yiV (fobsj-fmodelj) 



where fobs,i and f„mdeUi are the observed and corresponding 
model eigenfrequency of the ith mode, respectively, crohs,i is the 
uncertainty of the observed frequency, and A^ is the total num- 
ber of modes used for the fit. This approach provides a reliable 
and unbiased estimate of how well a set of observed frequen- 
cies coincides with a model eigenspectrum. A value of;^f^ < 1.0 
means that, on average, the model frequencies are within the un- 
certainty of the matched observed frequencies. For more details, 
see Guenther et al. (2005). However, if there is more than one 
model with;^'^ < 1.0, the set of observed frequencies (and un- 
certainties) is not sufficient to identify a unique model. In this 
case, additional stellar parameters (such as independently deter- 
mined temperature and/or luminosity) are necessary to narrow 
the parameter space of the model search. 

4.1. Mode Identifications 

Trying to match the observed oscillation spectrum to the mod- 
els is made significantly more difficult if we are not sure that 
all the observed frequencies are p-modes intrinsic to the star. 
With 59 frequencies identified as possible modes it is compu- 
tationally impractical to consider all possible combinations of 
p-modes versus spurious (i.e., not p-) modes. We therefor lim- 
ited our model search to the 21 most significant frequencies (sig 
> 10.0), as determined by SigSpec, in our list of 59 frequencies 
in the frequency range 20 to 95juHz. The searching code will try 
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Fig. 4. Same as Fig. |3] but for all photometric (circles) and ra- 
dial velocity (squares) frequencies. Symbol sizes are scaled with 
mode amplitudes: 40ppm to 125ppm and 0.9 to 3.4ms"' for the 
luminosity and radial velocity amplitudes, respectively. Arrows 
indicate RV modes which are presumably daily aliases. 



to find the best fit to as many frequencies as possible and this 
may, as is the case here, require that some modes be excluded 
from the fit. The modes excluded varies from model to model in 
the search. For the best fit model frequencies P23, P52, and P56 
are unmatched. We were unable to find any models whose ra- 
dial p-modes match more than 9 of the 21 observed frequencies. 
Only when we allowed the search to include nonradial modes up 
to and including / - 3 were we able to fit 18 (bold in Tab.[T]) of 
the 21 observed frequencies. The best model fit matches 9 of the 
frequencies to radial ;5-modes, 3 to / = 1, 5 to / = 2, and 1 to Z = 
3 modes with O-x^ - 0.22. The best fitting model has a mass of 
2.02 M0, an effective temperature of 4892K and a luminosity of 
60.1Lq, which places it within ±lcr of e Oph's position in the 
HR-diagram. In Table 2 we list the adiabatic ;?-mode frequencies 
of the best fitting model up to the acoustic cut-off frequency of 
~92fiHz. 

In Fig. [3] we shown, in an echelle diagram with folding fre- 
quency 5.2/iHz, the 21 most significant photometric frequencies, 
used in the search for the best fitting model, along with the p- 
mode frequencies of the best model fit. The filled circles corre- 
spond to the 18 photometric frequencies that are fit by the model. 
The p-mode frequencies of the best fitting model are shown as 
small grey filled symbols connected by line segments. In Fig.|4] 
we show both the MOST (circles) and the RV (squares) frequen- 
cies to provide a direct comparison between the two data sets. 
Needless to say, there are some observed frequencies that do not 
appear to fall on the Z-ridges of our best fit model. We argue be- 
low that some of these frequencies could still be intrinsic to the 
star 

The two frequencies P05 and P06 (6.44 and l.SA/Mz) that 
have the largest photometric amplitudes, but lie outside the 20 to 
95 yuHz range used to constrain the best fit model, do not match 
(within the uncertainties) any of the low n-valued /?-modes of the 
best fit model but are still close to the fundamental mode of the 
radial and 1=2 sequence. If these two frequencies are indeed low 
order p-modes than they imply that the deep interior structure of 
our models is sUghtly off. This is possible since the higher order 
modes mainly constrain the outer layers of the model. It is pos- 
sible that a different assumed initial helium abundance for the 
model, which would affect the frequencies of very lowest order 
modes, could yield a model that fits the lowest two frequencies 
along with the rest of the frequencies. We will address this and 
other modeling interpretations of the original 59 MOST frequen- 
cies in a separate paper 

Nine of the frequencies detected in the RV data match the 
p-mode frequencies of the best fitting model and are found in 
the list of photometrically observed frequencies, and 6 of the RV 
frequencies match the best fitting model and are not found in 
the list of photometrically observed frequencies. We note that 4 
RV frequencies would lie close to model frequencies if they are 
1 cycle/day aliases of the real frequencies. The presumed real 
frequencies are at the tips of the arrows in Fig.H] 

Although we are able to fit 18 of the 21 photometric fre- 
quencies with the largest SigSpec significances, only 63% ("*" 
in Tab.[lJ of the photometric frequencies (70% in the frequency 
region above 20juHz) and 76% of the RV frequencies can be at- 
tributed to pulsation modes in our best-fitting model. A few of 
the lowest frequencies may be associated with turbulence or ro- 
tational modulation in the envelope of the red giant or, as noted 
above, may not fit the model because the model does not have 
the correct helium abundance. We speculate, though, that fre- 
quencies that are not matched by the model could be the result of 
the damping and stochastical re-excitation of the modes. These 
modulations introduce spurious peaks around the true frequen- 
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cies of modes in Fourier space. We discuss this hypothesis in 
Section 4.4 below. 

We note that Hekker et al. (2006) attempted mode identifi- 
cations for the dominant frequency in the e Oph spectroscopic 
data (De Ridder et al. 2006). They identified the frequency R16 
at 58.2/iHz (having the largest RV amplitude) to be a mode of 
(I, m) - (2, 2). Our best fitting model matches this frequency with 
a mode of degree I - las well. 
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Fig. 5. Theoretical HR-Diagram showing the uncertainty box lo- 
cation of e Oph and a subset of the stellar model grid used for 
the pulsation analysis (light grey dots). The color scale gives the 
X^ values from the differences between the observed and model 
frequencies, where the scale is limited to values less than 1.0 
(for better visibility the color dark blue includes all models with 
X^ from 0.5 to 1.0). x^ < ^ indicated that the model frequencies 
are within 1 cr of the observed frequency uncertainty. The dark 
grey dots show models with a stellar radius of 1O.4+O.45R0 . 
The insert illustrates x^ as a function of model mass and age, 
respectively. 



the error box. We remind the reader that x^ is determined from 
the oscillation data only, that is, other observables are not used in 
any way to help determine the best model fit. The insets in Fig.|5] 
show;^'^ as a function of model mass and age. They illustrate the 
smooth distribution of x^ along the track of best fitting models 
(colored region) with constant large frequency separation. 

Richichi et al. (2005) list an angular diameter of 
2.94±0.08mas for e Oph in their catalogue of stellar radii ob- 
tained with the ESO VLTI. This angular diameter and the 
Hipparcos parallax of 30.34±0.79mas give a stellar radius of 
1O.4±O.45R0 . In Fig.|5] all models with a stellar radius within 
that range are marked by dark grey dots. The radius of the best- 
fitting model is about IO.8R0 , within the uncertainty of the in- 
terferometrically determined radius. 

4.3. Mode Amplitudes 

Using the effective temperature, mass and radius of our best- 
fitting model of e Oph (see Section 4.2), we obtain from the 
Kjeldsen & Bedding (1995) scaling relation a radial order of n 
= 9 at a frequency of ~59;uHz for the mode of maximum power 
The frequency with the largest RV amplitude is R14 at 58.2/iHz, 
corresponding to a mode with radial order « = 9 in our best fit- 
ting model and consistent with Kjeldsen and Bedding's scaling 
relation. 

Kjeldsen & Bedding (1995) provide a heuristically deter- 
mined relation between luminosity and RV amplitudes for solar- 
like oscillations scaled by the amplitude ratio of solar p-modes 
in light and velocity: 20.1 ppm per ms"' at a wavelength of 550 
nm. 

The average amplitude of the frequencies in the photome- 
try of e Oph matched to p-modes in the model (plus the two 
largest-amplitude peaks) is 85ppm. The corresponding value for 
the RV data is 1.64ms~^ This corresponds to an amplitude ra- 
tio of 51.8ppm per ms"'. If we perform a linear regression of 
the photometric and RV amplitudes for the subset of frequencies 
identified in both data sets, we obtain a slope of 41+7ppm per 
ms"' (we must neglect the frequency ~58//Hz with the largest 
RV amplitude to obtain this fit). This slope results in the follow- 
ing relation calibrated by the luminosity and temperature of e 
Oph and the central wavelength of the MOST passband: 



(SL/Lh 



Vosc/ms 



{AI550nm)(T,ffl5111Ky 



■26.4 ± 4.6ppm, 



4.2. Position ofe Opli in tine HR-Diagram 

A summary of fundamental parameters of the best fitting model 
and values found in the literature is given in Table 2. De Ridder 
et al. (2006) conducted an extensive literature search to locate 
e Oph reliably in the HR-Diagram. They adopted an effective 
temperature of 4887+lOOK and a luminosity of 59+5L0 . These 
values compare well with our best fit model which has a surface 
temperature of 4892K and a luminosity of 60.1 L0. 

In Fig.|5] we present a theoretical HR-Diagram showing a 
subset of the stellar model grid used for the pulsation analy- 
sis (small light grey points), the best model fits to the observed 
frequencies with x^ less than 1 .0 (color), and the error box for 
the luminosity and effective temperature of e Oph. The mod- 
els with the smallest deviation from the oscillation observations 
(0.28< x^ <0.3) are shown in red, and lie close to the centre of 



The coefficient (26.4±4.6 ppm per ms"') is close to that of the 
Kjeldsen & Bedding (1995) amplitude ratio for solar p-modes 
(20. 1 ppm per ms"'). However, we caution that the reduction of 
the e Oph spectroscopic time series was performed with respect 
to a nightly reference spectrum, which does suppress lower- 
frequency power and likely means a decrease of the measured 
RV amplitudes with increasing frequency relative to the intrinsic 
values. 

4.4. Mode Lifetimes 

Unlike p-modes in stars in the classical instability strip, which 
can coherently oscillate for millions of years, solar-type oscilla- 
tions driven stochastically by turbulence in the convective enve- 
lope have much shorter lifetimes. Pulsations of this nature can 
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be characterized by a superposition of randomly excited, intrin- 
sically damped, harmonic oscillations. The Fourier transform of 
a set of such incoherent signals produces peaks randomly dis- 
tributed under a Lorentzian profile with a width given by the 
damping rate (inversely proportional to the mode lifetime). 

In this work, we treat the individual frequencies in the data 
as independent radial and nonradial modes whose lifetimes are 
comparable to the duration of the observations. For solar-type 
oscillations, the scatter of observed frequencies about the true 
mode frequencies is a function of the mode lifetime (and the 
signal-to-noise ratio of the data). Kjeldsen et al. (2005) per- 
formed extensive simulations in their analysis of the RV oscil- 
lations of a Cen B to quantify the dependence of the frequency 
scatter on the lifetime and the S/N ratio. 

We apply a similar approach, assessing roughly 100,000 sim- 
ulated data sets using the time sampling of the MOST e Oph 
photometry. Each data set contains a damped and randomly re- 
excited oscillation at a given input frequency, generated by the 
method described by Chaplin et al. (1997), plus white noise of 
random amplitude. The generated mode lifetimes range from 0.5 
to 50 days. The average absolute deviation of the highest am- 
plitude peak (determined by SigSpec) from the input frequency, 
Av, as a function of the peak's SigSpec significance is presented 
in Fig. |6] for different mode lifetimes (solid lines). The left axis 
scale in this figure is specific to the MOST data run. From addi- 
tional simulations with data sets of different length T we found 
that the frequency deviation due to finite mode lifetime scales 
with the square root of T. The right-hand axis in Fig.|6]is cali- 
brated to be independent of T. We also present upper limits for 
the observational uncertainties of a stable mono-periodic sig- 
nal (horizontal dashed line scaled to the left axis) and a multi- 
periodic signal (diagonal dashed line, also, scaled to the left 
axis). (For more background details, see Kallinger et al. 2007b). 

As an example, consider some modes in our MOST obser- 
vations of £ Oph with SigSpec significance of 10 or more. If the 
mode lifetimes are 1 day on average, then according to our sim- 
ulation we can expect the observed frequencies of the modes to 
be on average within lyuHz of the intrinsic frequencies of the 
modes. If the modes have a lifetime of 20 days, then our sim- 
ulation indicates that the observed frequency of the modes will 
be on average within Q.l/Aiz of the intrinsic frequencies of the 
modes. 

We can invert this, by using the apparent scatter of observed 
frequencies to estimate the mode lifetimes. To determine the 
scatter in the observed frequencies we can either compare the 
frequencies of the modes observed in RV with the frequencies 
of the same modes observed in MOST photometry or we can 
compare the observed MOST frequencies with the nearest corre- 
sponding frequency in the best fit model. Unfortunately, the lat- 
ter method depends on the model fit itself. To avoid this we adopt 
the approach of Kjeldsen et al. (2005). We assume the modes 
to lie along ridges of common /-value and take the amount the 
modes deviate from the ridge as the scatter frequency (i.e., fre- 
quency difference). The frequency difference is taken to be the 
amount the mode's frequency deviates from a line drawn be- 
tween the mode below and the mode above the mode of com- 
mon /-value. Clearly, both of these methods depend on the cor- 
rect identification of the modes. We present the results of both 
methods. 

In Fig.|6]we plot, using light grey squares, the frequency dif- 
ferences (i.e., scatter) between the 12 modes common in the RV 
and MOST determined frequencies as a function of their SigSpec 
significance. The data points all (except one) lie below the 5d 
lifetime line of our simulation. Indeed, the average value of the 




30 40 50 60 



SigSpec significance 



Fig. 6. CaUbration of mode lifetimes for the observed oscilla- 
tion modes. As a result of numerical simulations of solar-type 
oscillations, solid lines indicate how much a detected frequency 
deviates on average from the input frequency as a function of 
the SigSpec significance for various mode lifetimes. The left axis 
gives the absolute deviation |Av| for the e Oph observing window 
and to be independent from the data set length T, the right axis 
gives the absolute deviation normalized to T'^^. The dotted lines 
correspond to the upper limit of the observational frequency un- 
certainty of a stable sinusoidal multi-periodic (horizontal line) 
and a mono-periodic signal (tilted line), respectively. Both are 
a function of T and are only valid for the left axis. Squares in- 
dicate the absolute deviations between frequencies determined 
from MOST and RV data for the same modes with an average 
value given by the large square symbol. Dots label the abso- 
lute deviation of all photometric mode frequencies from that ex- 
pected from the neighboring modes with equal degree /, using 
linear interpolation. Average values are shown as large dots for 
two ranges of SigSpec significance (below and above 10). 



frequency differences, shown as a black square in Fig.|6l corre- 
sponds to an average mode lifetime of ~12 days. 

Also in Fig.|6]we plot, using light grey circles, the frequency 
differences of the MOST photometric data assuming the intrin- 
sic mode frequencies do lie along smooth curves of common 
1-value. Although the scatter is greater than using the previous 
method, the averaged absolute deviations in the SigSpec signifi- 
cance range below and above 10 correspond to ~20 day and ~13 
day mode lifetimes, respectively (black circles in Fig.|6l). The 
light grey circles in the upper right corner of Fig.|6l correspond 
to low frequency modes, i.e., modes that naturally deviate the 
most from the regular spacing of the asymptotic formula. Both 
methods give similar average frequency deviations allocated to 
mode lifetimes of approximately 10 to 20 days. 

We also performed simulations to estimate the likelihood of 
a second ("spurious") peak appearing in Fourier space, given a 
mode with a lifetime 10-20 days and SigSpec significances com- 
parable to those observed in e Oph. The simulation assumed the 
time sample cadence of the MOST photometry and, separately, 
the time sample cadence of the RV data. We find that almost 
half of the frequencies will appear as "spurious" and apparently 
do not match the model, as we have found in Section 4. 1 when 
considering all 59 photometric frequencies. 
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5. Discussion & Conclusions 

By demonstrating the mutual consistency between our best fit 
models and the observations we conclude that observable nonra- 
dial modes do exist in giants. 

We have shown that the significant intrinsic frequencies 
present in the MOST photometry (Barban et al. 2007) and 
CORALIE + ELODIE RV measurements (De Ridder et al. 2006) 
of e Oph can be best matched to a model of radial and nonra- 
dial p-modes. The number of matched p-modes compared to un- 
matched modes as well as the scatter of the observed frequencies 
around the mode frequencies is consistent with mode lifetimes 
of 10-20 days. 

Like Barban et al. (2007), we find in the MOST data a clear 
well-populated sequence of radial modes with a fundamental 
spacing of ~ 5.8/iHz in the asymptotic regime and ~ 5.2/iHz 
in the range where most of the frequencies were observed. But 
for the 21 most significant frequencies in the photometric data 
above 20/iHz, we are only able to fit 9 of the frequencies to radial 
p-modes. Furthermore, we cannot find any model in our grid of 
more than 30,000 models that matches all 7 radial mode frequen- 
cies found by Barban et al. (2007) via Lorentzian profile fitting. 
We can, though, match 18 of the 21 most significant frequencies 
in the photometric data if we consider nonradial modes. 

The / - sequence is nearly fully populated and the / = 2 se- 
quence is also well populated, while the I - 1 and 3 sequences are 
less densely populated. This is not so surprising for I = 3 modes 
because their photometric and RV amplitudes in light integrated 
over the stellar disk are expected to be very small, even given the 
high precision of the data. But from the geometrical aspects of 
the spherical harmonics, I = 1 modes should have larger ampli- 
tudes than I - 2 modes for almost every orientation except nearly 
pulsational-equator-on. A possible explanation is suggested by 
Dziembowski et al. (2001). Their theoretical models show that 
1-2 nonradial modes are more frequently unstable than others, 
hence, / = 1 modes are less likely to be observable. In general, 
the distribution of modes according to /-value for e Oph resem- 
bles what Kallinger et al. (2007a) found for the red giant HD 
20884. 

The intrinsic nature of the most significant observed fre- 
quencies is also supported by the good agreement between 
the oscillation-only constrained best-fit model and the indepen- 
dently derived effective temperature, luminosity, and radius. Our 
best fit model has an effective temperature of 4892K and lu- 
minosity of 6O.1L0 compared to, e.g., De Ridder et al. (2006) 
non-asteroseismically derived temperature 4887±100K and lu- 
minosity 59±5L0. Furthermore our best model fit has a radius of 
1O.82R0 that agrees with the radius 1O.4±O.45R0 based on ESO 
VLTI angular diameter determination of Richichi et al. (2005) 
and the Hipparcos parallax. The mass and age of e Oph as con- 
strained by the oscillations is 2.O2M0 and 0.77Gyr. 
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Table 1. Oscillation frequencies (f), amplitudes (A) and SigSpec 
significances (sig) for e Oph. Frequencies used to constrain the 
fit are marked in bold face. P23, P52 and P56 are also used in the 
model search procedure but do not match any mode of the best 
fitting model. "+" or "-" indicate RV frequencies which are pre- 
sumably + 1 d"' or -1 d ' aliases. Asterisk symbols indicate fre- 
quencies matched by modes of the best fitting model (see Tab.|2]i. 
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Table 2. Adiabatic ;?-mode frequencies of the best fitting model 
up to the acoustic cut-off frequency in units of yuHz. For nonra- 
dial modes Up is the number of radial nodes in the /9-mode part 
of the modes. Identifier of observed modes matching a model 
frequency are listed in brackets. Bold faced modes indicate the 
modes matched by observed frequencies used for model fitting. 
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6.797 (™^) 
13.856 P 
20.119 P 

26.062 
31.549 p 

36.562 ("3) 
41.892 ("') 
47.002 (''!«) 
52.243 P 
57.693 P 
63.074 P 
68.508 C^J) 
74.035 
79.559 
85.086 P 
90.625 



8.416 P 
15.561 P 
21.807 {^^) 
27.779 (™) 
33.089 (™) 
38.235 ("') 
43.585 P 
48.694 P 
54.078 P 
59.539 
64.940 
70.434 
75.985 
81.516 
87.061 
92.597 



Table 3. Fundamental stellar parameters for e Oph as found in 
the literature and for the best fitting model. 



Literature 



best fitting 
model 



Effective temperature ... [K] 

Luminosity ... [Lq] 

Radius ... [Rq] 

Mass ... [Mq] 

logg... [gems"-] 

Age ... [Gyr] 

Mixing length parameter ... [Hp] 

Metallicity ... [Z] 



4877±100' 


4892 


59±5' 


60.13 


10.4±0.452 


10.82 




2.02 


2.48±0.363 


2.674 




0.770 




1.74 


0.01* 


0.01 



'DeRidderetal. (2006) 

^Richichi et al. (2005) and Hipparcos 

^Allende & Lambert (1999) 

"[Fe/H] = -0.25 from Cayrel de Strobel et al. (2001) 



